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O \ Abstract: We address the problem of decomposing graphs in perturbative QCD 

i-^ \ into terms associated with particular regions. Motivated by asking how to incor- 

I \ porate next-to-leading order (NLO) QCD corrections in parton-shower algorithms, 

D \ we require that: (a) The integrand for the hard part is to be integrable even if 

'T^ \ the corrections are applied to a process that is not infrared and collinear safe, (b) 

.^ \ The splitting between the terms should be defined gauge-invariantly. (c) The de- 

r> \ pendence on cut-offs should obey homogeneous evolution equations. In the context 

c^ \ of one-gluon-emission graphs for deep inelastic scattering, we explain a subtractive 

technique that is based on gauge-invariant Wilson-line operators. Appropriate or- 
ganization of subtractions involving the soft region allows a connection to previous 
work where evolution equations with respect to the directions of the Wilson lines 
have been derived. 
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1. Introduction 

Current parton-shower Monte Carlo event generators are essentially leading-order 
(LO) QCD tools. For precision phenomenology at present and future high-energy 
colliders, it is valuable to be able to go beyond this level of approximation [1]. Al- 
though there are a number of treatments of various next-to-leading-order (NLO) 
effects in event generators (see [2, 3, 4, 5] and references in [1]), there is as yet no 
method for going beyond the leading approximation systematically. This implies 
that event generators are not able to incorporate fully the known NLO (and NNLO) 
calculations of hard scattering cross sections. 

Recently, systematic subtractive procedures have been proposed [6, 7, 8] to cor- 
rect this situation. An important step in the implementation of this program is to 
show how to decompose Feynman graphs into sums of terms over different regions; 
the terms are to be arranged to correspond to factors in a factorization theorem that 
is in a form suitable for the Monte Carlo application. In Ref. [7] an implementation 
was presented for a simple, but phenomenologically important case: the photon- 
gluon fusion process in leptoproduction. This case was simple because soft gluons do 
not enter at the leading power, so that leading regions do not overlap. The purpose 
of the present work is to show how to extend the method to decompose graphs with 
soft gluons and hence with overlapping leading regions. We will treat the simplest 
such case in leptoproduction, i.e., the photon-quark process to NLO. 

A Monte Carlo event generator treats observables that are not infrared and 
collinear safe. So the kind of subtraction method used in Refs. [9] and [10] is not 
directly applicable, since it relies on the calculated observable being infrared and 
collinear safe. In particular we cannot use a cancellation of soft gluon contributions 



between real and virtual graphs. More general factorization theorems are to be used 
in which the soft contributions factor instead of canceling [11], and all the factors 
are defined gauge-invariantly in terms of Wilson-line operators [12]. We will require 
that the following properties be satisfied: 

(a) The integrand for the hard scattering coefficient is to be an integrable function, 
and not merely a distribution. Thus the hard scattering coefficient is usable 
even when one is calculating an infrared and collinear unsafe observable. See 
Eq. (3.9) for an example. 

(b) The terms in the expansion of each Feynman graph should arise from matrix 
elements of gauge-invariant operators. 

(c) In particular, the necessary cut-offs on rapidity integrations are to be defined 
gauge-invariantly. This involves the use of Wilson-line operators along non- 
lightlike directions. 

(d) The evolution equations [12, 13] with respect to these cut-offs should be sim- 
ple, in the sense that there should be no power-law remainder terms. That 
is, the equations are homogeneous, like the renormalization-group equations, 
rather than the Callan-Symanzik equations. This is achieved by the technique 
proposed in Ref. [14]. 

In Sec. 2 we briefly describe the framework provided by the subtraction method 
for Monte Carlo event generators. In Sec. 3 we treat gauge-invariant subtractions 
and construct the decomposition of the partonic cross section. In Sec. 4 we comment 
on graph-by-graph subtractions. In Sec. 5 we give conclusions. 

2. Subtraction scheme 

To put our result in the Monte Carlo context, we schematically represent the cross 
section in an event generator as 

a[W]= Y. W{X)V^®H. (2.1) 

final states X 

Here H^ is a weight function that specifies the definition of the particular cross section 
under consideration. The symbol PS denotes the parton shower and the symbol ® 
denotes its action on the initial and final partons in the hard scattering, whose cross 
section is denoted by H . In a standard Monte Carlo, the hard scattering is taken to 
the leading order, iJ^^^), and PS denotes showering from the partons in H^^'~'\ In 
an NLO Monte Carlo, the cross section involves a structure of the form 



PS 



i7(LO) + ^^ fjj(^LO) _ pg^(^) ^ ^(LO) _ pg^^i) ^ ^(LO)\l ^ ^2.2) 



Here the first term in the square brackets is the LO hard-scattering function, and the 
second term is the subtracted NLO hard-scattering function, if (^^*^) is the resuh of 
computing the partonic cross section from the NLO graphs, while PS/(1) and PSir(l) 
are the order a^ approximations to the initial-state and final-state showering. The 
subtraction terms avoid double counting of events already included by showering 
from ff(LO) 

We will construct a decomposition of H^^^'^^ into a sum of pieces, one for each 
of the leading regions R, 

j^(NLO) ^ J2 Ah{R) + nonleading power, (2.3) 

regions R 

that holds uniformly over the whole of the phase space. Each of the pieces in (2.3) 
will contain counterterms that prevent double counting and provide the suppression 
for going outside the region in which that particular piece was originally supposed 
to give a good approximation to the matrix element. This is to be contrasted with 
approaches based on splitting the phase space in different domains and using different 
approximations to the matrix element in these different domains (see, e.g., Refs. [2, 
3]). 

Once we have established a result of the type (2.3), the term corresponding to 
the ultraviolet region will give precisely the subtracted hard-scattering function to be 
used in Eq. (2.2). The collinear terms will correspond to factors in the cross section 
that are associated with showering; they therefore imply the evolution kernels to be 
used in the showering. The soft term would correspond to a new element in the 
Monte Carlo, but, at the order to which we work, we will show that the soft term 
can be eliminated by a suitable choice of the cut-offs. This is a result analogous to 
one in [9]. However the results of [11] indicate that this is unlikely to be an all-order 
result. 

Our decomposition of graphs of order Ug entails a specific definition of the 
collinear factors, which will not necessarily coincide with the definitions used in any 
current event generator. However, we will not address in this paper the issue of how 
to obtain a showering algorithm that corresponds to our collinear terms. A correct 
solution of this problem will encompass known results about coherent emission of 
gluons and angular ordering of gluon emission [15]. 

3. Decomposition of partonic cross section 

Let us consider deep inelastic scattering 7* -|- P ^ X, in which we consider a generic 
cross section or observable associated with the reaction 7*(g) -|- q{p) -^ g{k) + q{k') 
— see Fig. 1. We denote this by S [(/?], where (/? is a weight function that is the 
product of the weight function W{X) concerned with the final state X in Eq. (2.1) 
and the factors in the cross section associated with the showering, including the 



parton density. Thus if contains all the infrared sensitive and nonperturbative parts 
of the observable. 
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Figure 1: DIS with quark-induced hard scattering. 



We work in the 7* + hadron reference frame, and we use light-front coordinates 



{v^ ^v , vt) with v^ 



± v'^)/ \/2. The hadron and photon momenta are 



pi^ = (P+^ m^/2P+, Ot) and q^' = {-xP+, Q'^/2xP+, Or). Then we parameterize the 
gluon momentum k as 



*'=H*"^2S^ 



\^t\ 



(3.1) 



where is a unit transverse vector at azimuthal angle 0. 

For our calculation, the external partons are on-shell, and the incoming quark p 
has zero transverse momentum, so that S can be written as follows: 



S[v.] 



00 /-oo /•27r (j^A) 
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Here, J is the Jacobian factor 



J(x, a, P) 



1 



1 



167r2 l + a-P 
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(3.2) 



(3.3) 



and Ai is the next-to-leading-order matrix element for 7*g obtained by contracting 
the photon Lorentz indices with the projector corresponding to the structure function 
F2 [16] 



M = 4e'^ gl Cf M{a, P) , M(a, P) = {I - PY 



l + (l+a-pY 
a P {1 + a- p) 



2 + 6 



[i-PY 



l + a- p 
(3.4) 

The physical region for a,P is the interior of the triangle in Fig. 2. 

Standard arguments [17] determine the infrared sensitive regions contributing 

to the leading power behavior of S[(p], which are located on Fig. 2 as follows: The 

region in which the gluon is collinear to the initial state is a neighborhood of the axis 




(l-x)/x 

Figure 2: The phase space of Eq. (3.2) in the a, f3 plane. 

(3 = 0, the region in which the gluon is colhnear to the final state is a neighborhood 
of the axis a = 0, and the soft region is a neighborhood of the origin a = 0, j3 = 0. 
The truly hard region lies away from the a = and (3 = axes. 

To obtain a decomposition for S of the type of Eq. (2.3), we now employ the 
technique of Ref. [14]. This generalizes the R-operation techniques of renormal- 
ization. (See Ref. [18] for a related approach.) To ensure that the procedure is 
gauge-invariant, each of the terms in the right hand side of Eq. (2.3) is constructed 
from matrix elements involving Wilson line operators, 

Vi{n) = V exp iig dyn-A{yn)] , Vp{n) = V exp iig dyn-A{yn)], 

(3.5) 
with suitable directions n for the lines. Evolution equations in n enable one to 
connect the results corresponding to different directions [12, 13, 14]. We define light- 
like vectors p = (1,0, Ot), p' = (0,1,07). We will also use non-lightlike vectors 
u = {u~^,u~,Ot), u' = {u'~^,u'~ ,0t), all of whose components are positive. It is 
convenient to define rj = (2x^P+ /Q'^)u~ /u^ , and rj' = {Q"^ /2x'^P'^ )u'^ /u'" . 

As in [14], we start with the smallest region, the soft region a,/? — > 0, and 
determine the corresponding contribution to the matrix element (3.4): 

Ms{a, /5) = 4 - ( Jm^ R - (R^ V (3.6) 

ap [a + r]'(i) p a [p + rja) 

Observe that the first term in the right hand side of this formula is just obtained 
by taking the soft approximation to Eq. (3.4). It can be thought of as the one- 
loop contribution to the square of a vacuum-to-gluon matrix element of a product 
of eikonal Wilson lines taken along lightlike directions p, p' [14]. This first term 
reproduces the behavior of the matrix element M when a and (5 simultaneously 
approach zero. But there are also logarithms in its integral associated with the 
collinear regions where aj (5 or /3/q; go to zero. The subtractions provided by the other 
two terms conveniently cancel these regions. They can be derived from operators 
analogous to those for the first term, except for replacing one of the lightlike eikonal 
lines by a line along a non-lightlike direction. In particular, the second term subtracts 



the divergence from a region collinear to the initial state, i.e., from the region (3 /a -^ 
0. At the same time, the non-hghthke vector u' in this term provides a cut-off on 
the region of small a. Similarly, the third term in Eq. (3.6) subtracts the divergence 
from the region collinear to the final state, i.e., from the region a/ {3 -^ 0. 

Next we construct terms for the collinear regions. By applying a treatment [14] 
analogous to that for the soft region, we arrive at 

,w m 1 l + (l + a)^ 2 2 




p a (l + a) aj3 {a + r]'p) p 
for the region collinear to the initial state, and 

M (r. R\ 1 (l-/3) + (l-/3f 2 2 

a p ap a [p + rja) 

= -(-4 + 3/3-/5^ + -^) (3.8) 

a p + rja 

for the region collinear to the final state. The first term of the expression in the 
middle in each of these equations is the unsubtracted collinear approximation to the 
original matrix element (i.e., the f3 -^ ot a ^ limit of M). We will comment 
below on the subtraction terms. 

The fully subtracted matrix element, associated with the hard region, is then 
given by 

Mnia, f3) = M - Ms - Mj - Mp 



a (3 {1 + a- (3) I + a - (3 

_ 1 l + {l + ay _ 1 1 + (1-/J)^ (1 _ ^) + J_ 
P a {1 + a) a /3 a P 

3, ^ I 6(^-^)' (39) 

^ {I + a){l + a - p) {I + a ~ py ^ ' ' 

This matrix element is finite in all of the infrared regions. It can be safely integrated 
down to a = or /3 = 0. Moreover, it is independent of the choice of the non-lightlike 
directions -u, u': the dependence on 77, 77' has canceled in Eq. (3.9). 

Eqs. (3.6)-(3.9) achieve a decomposition of the type (2.3) for the original matrix 
element. There is one term for each of the leading regions — in particular, a soft 
term. We now ask: can we reorganize it in a way that is suited for a parton-shower 
algorithm, such as, e.g., the algorithm [19] used in the event generators [20]? 

The soft term can be eliminated by choosing the vectors u and u' so that 7777' = 1; 
then there are only collinear terms, as is appropriate to match the structure of 
the parton-shower Monte Carlo algorithm, which has only initial-state or final-state 
branchings. The symmetric choice 77 = 77' = 1 gives 

M<MC)(„_^),l l + (l+af _l 2 

^ ^ ' pa (1 + a) a a + p ^ ' 



Afr'(a./?) = iii^^(l-/3)-i-^. (3.11) 

a p p a + p 

All of the infrared contributions are now associated with configurations that are ei- 
ther collinear to the initial state or collinear to the final state. We have inserted 
superscripts in the left hand sides of Eqs. (3.10) and (3.11) to indicate that this par- 
ticular choice of the non-lightlike directions gives rise to a structure that corresponds 
to that of Monte Carlo algorithms. 

Eq. (3.9) determines the form of the counterterms to be used in the matrix 
element; Eqs. (3.10) and (3.11) determine the form of the counterterms to be used in 
the showering. The subtracted NLO matrix element (3.9) only receives contributions 
from the truly ultraviolet region. As for the modification to the showering, consider 
Eq. (3.10). M) is associated with the first branching from the initial state. As 
noted below Eq. (3.8), the first term in the right hand side is obtained from taking 
the collinear approximation /? ^ to the original matrix element. The coefficient of 
1//9 is the standard quark — » quark splitting kernel: 



a{l + a) ""Vl + a 

This first term corresponds to the standard form of the showering. It gives a good 
approximation in the initial-state collinear region, i.e., /5 ~ 0. The second term in 
the right hand side of Eq. (3.10), on the other hand, is non-standard. It effectively 
provides a cut-off when a ^ 0. Note that the second term is suppressed in the 
collinear region /5 ^ at fixed a. That is, the modified showering coincides with 
the usual one in the collinear limit and differs from it away from the collinear limit. 
Analogous remarks can be made based on the formula (3.11) associated with the 
final state. 

Observe that if one regulated the a — > behavior of the first term in Eq. (3.10) 
by subtracting its a — > limit, given by 2/{aP), this would bring about an extra 
/9 — s> singularity. This would not be suited for our application in a Monte Carlo 
algorithm. In contrast, the second term in Eq. (3.10) represents precisely what is, 
from our point of view, a better choice of a counterterm: it subtracts the a ^ 
singularity without introducing any extra singular behavior at (3^0. 

Note also that this counterterm cuts off the integration over the region of small 
a at a value of order /3: the leading behavior of Mj for small a is of the type 

o 

-|- regular terms in a. (3.13) 



(3 (a + (3) 



Then we see that the procedure based on gauge-invariant subtractions that we have 
just applied, compared to the cut-off method, tells us precisely where the cut-off is 
to be placed. The position of the cut-off on a turns out to be /3-dependent. In more 
physical terms, this indicates that the cut-off to be applied in the initial-state shower 



and the cut-off to be applied in tlie final-state shower are not to be set independently, 
but they are related. 

4. Graph-by-graph subtractions 

In our calculational example, the form of the four terms was determined by requiring 
them to be obtainable from matrix elements of gauge-invariant operators. However, 
the procedure can be applied graph-by-graph, as we will now explain. The basis of 
this procedure is the derivation of factorization for soft factors given in Refs. [11, 12]. 
Consider the connection of a gluon to a subgraph that consists of lines that are 
all collinear in the + direction. This factor we denote by J'^{k), where n is the 
Lorentz index that couples to the gluon. In a region where the gluon is either soft or 
collinear in the opposite direction, it is a good approximation to replace J'^ by its + 
component, and to replace the momentum /c in J by its — component. Multiplying 
and dividing by A;~ gives 

= ^'^i-[fc-J+(0,A;-,OT)]. (4.1) 

Since the last factor is of the form of a gluon Green function contracted with the 
gluon's momentum, a Ward identity can be applied. Then after a sum over all 
relevant graphs we obtain factors corresponding to matrix elements of operators that 
include path-ordered exponentials [12]. 

When this replacement is made, together with the corresponding replacement for 
gluonic connections to the opposite "jet" factor, some of the terms have new diver- 
gences when the gluon's rapidity goes to +oo or — oo, as explained earlier and in Ref. 
[14]. These divergences are cancelled by further subtractions that are constructed by 
manipulations on the 1/k^ and 1/k'^ factors, to give a result like Eq. (3.6). 

5. Conclusions 

In conclusion, the above calculation tells us how to organize the subtractions in situa- 
tions in which soft gluons are present and both initial-state and final-state branchings 
contribute. This is one of the issues that have to be dealt with to construct NLO 
Monte Carlo event generators. The subtractive technique [14] used for this calcula- 
tion is based on graph-by-graph R-operation methods and enforces gauge invariance 
by relating the counterterms to matrix elements of Wilson-line operators. The cal- 
culations of this paper do not address the issue of how to generate the whole shower 
corresponding to the subtracted collinear terms. The investigation of this is left to 
future work. 
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